Detecting barrier to cross-jet Lagrangian transport and its destruction in a 

meandering flow 
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Cross-jet transport of passive scalars in a kinematic model of the meandering laminar two- 
dimensional incompressible flow which is known to produce chaotic mixing is studied. We develop a 
method for detecting barriers to cross-jet transport in the phase space which is a physical space for 
our model. Using tools from theory of nontwist maps, we construct a central invariant curve and 
compute its characteristics that may serve good indicators of the existence of a central transport 
barrier, its strength, and topology. Computing fractal dimension, length, and winding number of 
that curve in the parameter space, we study in detail change of its geometry and its destruction that 
are caused by local bifurcations and a global bifurcation known as reconnection of separatrices of res- 
onances. Scenarios of reconnection are different for odd and even resonances. The central invariant 
curves with rational and irrational (noble) values of winding numbers are arranged into hierarchical 
series which are described in terms of continued fractions. Destruction of central transport barrier 
is illustrated for two ways in the parameter space: when moving along resonant bifurcation curves 
with rational values of the winding number and along curves with noble (irrational) values. 

PACS numbers: 05.45.-a,05.60.Cd,47.52.+j 



I. INTRODUCTION 

A meandering jet is a fundamental structure in labo- 
ratory and geophysical fluid flows. Strong oceanic and 
atmospheric jet currents separate water and air masses 
with distinct physical properties. For example, the Gulf 
Stream separates the colder and fresher slope ocean wa- 
ters from the salty and warmer Sargasso sea ones. Re- 
cently, there has been much interest in applying ideas and 
methods from dynamical systems theory to study mixing 
and transport in meandering jets. In steady horizontal 
velocity flelds, water (air) parcels move along streamlines 
in a regular way. When the velocity fleld changes in time, 
the motion becomes much more complicated even if the 
change is periodic. The phenomenon of chaotic advection 
of passive particles in (quasi) periodically-disturbed fluid 
flows has been studied theoretically and experimentally 

In the context of dynamical systems theory, chaotic ad- 
vection is Hamiltonian chaos in two-dimensional incom- 
pressible flows (for a review of Hamiltonian chaos see, 
for example, [1,0, 0, Q)- The coordinates x and ?/ of a 
passive particle on the horizontal plane satisfy to simple 
Lagrangian equations of motion 
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where u and v are zonal and meridional velocities of the 
particle, and the stream function ^' plays the role of a 
Hamiltonian. The phase space of the dynamical system 
(HJ) with one and half degrees of freedom is a configuration 
space for passive particles. 

A number of simple kinematic and dynamically consis- 
tent model stream functions have been proposed to study 
large-scale chaotic mixing and transport in geophysical 



meandering jet fiows [HI, [13, [H, 0, E [3 - De- 

terministic models do not pretend to quantify transport 
fluxes in real oceanic and atmospheric currents but they 
are useful to reveal large-scale space-time structures that 
specify qualitatively mixing and transport of water and 
air masses. Whether or not the jet provides an effective 
barrier to meridional or cross-jet transport, under which 
conditions the barrier becomes permeable and to which 
extent, these are crucial questions in physical oceanogra- 
phy and physics of the atmosphere. The problem must be 
treated from different points of view. In the straightfor- 
ward numerical approach based on full-physics nonlinear 
models, the velocity field is generated as an outcome of a 
basin circulation model and a flux across the jet (if any) 
can be estimated integrating a large number of tracers. 
The kinematic and linear dynamically consistent models 
are less realistic, but they allows to identify and ana- 
lyze different factors which could enhance or suppress 
the cross-barrier transport. 

As to meandering currents, both the approaches have 
been applied to study cross-jet Lagrangian transport. A 
simple kinematic model with the basic streamfunction 
in the form ^ has been shown to reproduce some fea- 
tures of the large-scale Lagrangian dynamics of the Gulf 
Stream water masses The phase portrait of Eqs. ([T]) 
with the meandering Bickley jet ^ is plotted in Fig.[T](a) 
in the frame moving with the meander phase velocity. 
Time dependence of the meander amplitude or the in- 
troduction of a secondary meander, superimposed on the 
basic flow, may break the boundaries between distinct 
regions in Fig, [l] (a) producin g ch aotic mixing and trans- 
port between them [a, [13, [13, lli, [H, [H • The numerical 
calculations, based on computing the Melnikov function 
[20| . have shown that transport across the jet was much 
weaker than that between the jet (J), the circulation cells 
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(C), and the peripheral currents (P) in Fig. [T](a), i. e. 
the perturbation mixes the water along each side of the 
jet more efficiently than across the jet core An at- 
tempt to analytically predict the parameter values for the 
destruction of the transport barrier was made in Ref . [l^ 
using the heuristic Chirikov criterion for overlapping res- 
onances A technique, based on computing the finite- 
scale Lyapunov exponent, as a function of initial position 
of tracers, has been found useful in Ref. [l^l to detect the 
presence of cross-jet barriers in the kinematic model |[2j- 
An analysis of cross-jet transport, based on lobe dynam- 
ics, has been applied in [23| to describe how particles can 
cross the jet from the north to the south and vice versa. 

The study of cross-jet transport has been motivated 
also by a series of laboratory experiments [13, HI, 
on Rossby waves propagating along an azimuthal jet in 
a rapidly rotating tank. This flow can be modeled in 
the Hnear approximation of the corresponding fluid equa- 
tions [131 by a stream function which is a superposition 
of a Bickley jet and two neutral modes (Rossby waves). 
The destruction of a barrier to cross-jet transport has 
been studied analytically by using the Chirikov criterion 
in the pendulum approximation and numerically by us- 
ing Poincare sections j23|. It was shown that one needs 
very large values of the perturbation amplitudes to break 
the barrier. The analytic model proposed in [131 has 
been used recently to study Lagrangian dynamics of at- 
mospheric zonal jets and the permeability of the strato- 
spheric polar vortex. Poincare sections and flnite-time 
Lyapunov exponents revealed a robust transport barrier 
which can be broken either due to large perturbation am- 
plitudes of the Rossby waves or as a result of an increase 
of their phase velocities A comparison of properties 
of cross-jet transport in ad hoc kinematic and dynami- 
cally consistent models of atmospheric zonal jets has been 
done recently in Ref. (29| . 

Being motivated by Lagrangian observations of the 
oceanic currents, cross-jet transport and mixing have 
been studied in numerical models of meandering jets 
[ill . It has been shown both in barotropic 

and baroclinic nonlinear numerical models, where the 
meander amplitude can not be made arbitrary large, 
that cross-jet chaotic transport, resulting from the me- 
andering motions, are maximized at a subsurface level. 
Since the undisturbed velocity is weaker at deeper lev- 
els, the corresponding separatrices are closer to the jet 
core. Therefore, separatrix reconnection should occur 
below some critical depth and transport across the jet 
should be facilitated. 

Independent on the work on cross-jet transport in 
the geophysical community, there have been a number 
of theoretical and numerical investigations of chaotic 
transport in so-called area-preserving nontwist maps 
[13, nil, M, m, 113, m, m, IB, we mention spe- 

cially the early study of different reconnection scenario 
[31, 3§| and the flrst systematic study of cross-jet trans- 
port in nontwist maps [13, [13, [1H| . These maps locally 
violate the twist condition, a map analogue of the non- 



degeneracy condition for Hamiltonian systems. Nontwist 
maps are of interest because many important mathemat- 
ical results, including KAM and Aubry-Mather theory, 
depend on the twist condition. Apart from their mathe- 
matical importance, nontwist maps are of a physical in- 
terest because they are able to model transition to global 
chaos, the term meaning in the mathematical community 
a cross-jet transport. Nontwist maps allow to study dif- 
ferent scenarios for this transition: reconnection of sep- 
aratrices, meandering and breakup of invariant tori, and 
others. 

The onset of global chaos in oscillatory Hamiltonian 
systems, where the eigenfrequency possesses a local ex- 
tremum as a function of energy, has been studied analyt- 
ically and numerically in Refs. [H, |43|- In such systems 
with two or more separatrices, global chaos may occur at 
unusually small magnitudes of perturbation due to over- 
lap in the phase space between resonances of the same 
order and their overlap in energy with chaotic layers of 
the corresponding unperturbed separatrices. 

In the present paper we develop a method for detect- 
ing a barrier to cross-jet Lagrangian transport (or global 
chaos in a more general context), apply it to the kine- 
matic model of a meandering jet flow, study changes 
in its topology under varying the perturbation param- 
eters, and scenarios of its destruction. In Sec. II we 
briefly introduce a model streamfunction which is known 
to produce chaotic advection [1, [H, [IS] and compute the 
amplitude-frequency s-v diagram demonstrating the pa- 
rameter range for which cross-jet transport exists. Based 
on the symmetry of the flow, we propose in Sec. Ilia a 
numerical method to identify a central invariant curve 
(CIC) which is a diagnostic means to detect the process 
of destruction of a central transport barrier (CTB). The 
CIC is constructed by successive iterations of so-called in- 
dicator points jl3| . Computing the fractal dimension of a 
set of iterations of those points at different values of the 
parameters, we identify whether the CIC and CTB are 
broken or not. In Sec. Illb we study possible geometries 
of the CIC that may change dramatically with varying e 
and ly. Before the total destruction, the CIC experiences 
a number of local bifurcations becoming a complicated 
meandering curve whose properties can be specifled by 
its length and the winding number w. A structure of 
the set of CICs is revealed in a continued fraction rep- 
resentation of their winding numbers. The CICs with 
rational w are arranged in hierarchical series connected 
with the corresponding resonances. Whereas, the CICs 
with noble numbers form their own series. Destruction 
of CTB is studied in Sec. IV for two ways in the pa- 
rameter space. When moving along a so-called resonant 
bifurcation curve with a rational value of w, one specifles 
the values of e and v for which the CIC is broken but 
CTB remains. In contrary to that, when moving along 
any curve with noble value of w, a CIC exists providing 
CTB. The process of CTB destruction in both the cases 
is illustrated in Sec. IV. 
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II. THE AMPLITUDE FREQUENCY 
DIAGRAM FOR CROSS-JET TRANSPORT IN 
THE MODEL FLOW 

We take the Bickley jet with a running wave imposed 
as a kinematic model of a meandering shear flow in the 
ocean. The respective normalized stream function in the 
frame moving with the phase velocity of the meander has 
the following form 



tanh 



y — A cos X 



Cy, 



(2) 



where the jet's width Ljet, meander's amplitude A and 
its phase velocity C are the control parameters. The 
phase portrait of the advection equations |(T]) with the 
streamfunction shown in Fig.[T](a), consists of three 
different regions: the central eastward jet J, chains of 
the northern and southern circulation cells C and the 
peripheral westward currents P. The flow is steady in the 
moving frame of reference, and passive particles follow 
the streamlines. In Fig. [T] (b) we plot a frequency map 
f{xo,yo) that shows by nuances of the grey color the 
value of the frequency / of particles with initial positions 
(xo,yo) in the unperturbed system. The maximal value 
of the frequency, /max = 1.278, have the particles moving 
in the central jet. 

As a perturbation, we take the simple periodic modu- 
lation of the meander's amplitude 



A — Ao + ecosvt. 



(3) 



Under the perturbation, the separatrices, connecting sad- 
dle points, are destroyed and transformed into stochastic 
layers. The strength of chaos depends strongly on both 
the perturbation parameters, the perturbation amplitude 
e and frequency v. In the model used the normalized con- 
trol parameters are connected with the dimensional ones 
as follows [13]: A = afc, C = c/uj„Xk, and Ljet = Afc, 
where a, k, and c are amplitude, wave number, and phase 
velocity of a meander, respectively, A and are char- 
acteristic width and maximal zonal velocity in the jet 
on the surface. All these parameters change in a wide 
range in the Gulf Stream [i, Oj, |li|: A ~ 40 100 km, 
a ~ 50^60 km, 27r/fc ~ 200H-400km, c ~ 0.1^0.5 m/sec. 
Mm ~ 1-^1.5 m/sec. So, we get Ljet — 0.1-^3, A ~ 0.7-^2, 
and C ~ 0.02 0.3. Being motivated by mixing and 
transport in the Gulf Stream, we took the following nor- 
malized values of the control parameters that will be used 
in all our numerical experiments: Aq — 0.785, C — 0.1168 
and Ljet = 0.628. 

The equations of motion (U) with the stream function 
([2]) and the perturbation |(3]) have the symmetry 



S 



y 



-y 



and the time reversal symmetry 

y' = y- 



^0 



(4) 



(5) 



The symmetries ([4]) and ([5]) are involutions, i. e. — 1 
and Iq = 1. Due to the symmetry S, motion can be 
considered on the cylinder with < a; < 27r. The part 
of the phase space with 2?™ < x < 27r(n + 1)^ n — 
0, ±1, . . . , is called a frame. It should be stressed that 
the phase space in two-dimensional incompressible flows 
is a configuration space for advected particles. 




FIG. 1: (a) Phase portrait of the unperturbed flow in the 
frame moving with the meander's phase velocity. Streamlines 
in the circulation cells (C), jet (J), and peripheral currents 
(P) are shown, (b) Frequency map represents by color values 
of the frequency / of particles with initial positions {xo,yo) 
advected by the unperturbed flow. 

The following numerical procedure has been applied 
to estabhsh the fact of cross-jet transport in the kine- 
matic model of the meandering jet current. The ad- 
vection equations |(T|) for given values of the perturba- 
tion amplitude e and frequency v and with twenty par- 
ticles, released nearby the northern saddle point, have 
been integrated up to the time instant when one of the 
particles was detected to cross the straight line y — ys 
passing through the southern saddle. If after the time 
Lmax = 1000 x2n/v none of the particles crosses the Hne 
y = ys, we assume that for given values of the parame- 
ters there is no cross-jet transport. The e-v diagram in 
Fig. [2] shows the values of the parameters for which cross- 
jet chaotic transport exists (white-color zones). There 
are a number of the frequency values for which transport 
occurs at surprisingly small values of the perturbation 
amplitude e. The absolute minimal value of the pertur- 
bation amplitude, at which the cross-jet chaotic transport 
occurs, Eniin = 0.0218 « Aq/36, corresponds to the fre- 
quency V = 1.165 which is close to the natural frequencies 
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of the particles moving in the central jet. 



V 
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e 

FIG. 2: Amplitude-frequency e-u diagram showing the pa- 
rameter values for which cross-jet chaotic transport exists 
(white zones) or not (black zones). 



III. TOPOLOGY OF A BARRIER TO 
CROSS-JET TRANSPORT AND CENTRAL 
INVARIANT CURVE 

The amplitude-frequency diagram is useful to detect 
cross-jet transport but its computation is a time con- 
suming procedure. Moreover, it says nothing about the 
properties of barrier to transport and mechanism of its 
destruction. A fractal-like boundary between the colors 
in Fig. [2] reflects an intermittency in appearance and de- 
struction of the cross-jet barrier when varying e and v. 
Further insight into topology of the barrier could be ob- 
tained if one would be able to find an indicator of cross- 
jet transport, i. e. an object in the phase space whose 
form contains an information about permeability of the 
barrier. 



A. Detecting the central invariant curve 

First of all, we need to give definitions of some basic 
structures specifying a cross-jet barrier and its destruc- 
tion. The central transport harrier (CTB) is defined as 
a strip between the southern and northern unperturbed 
separatrices confined by marginal northern and southern 
baUistic trajectories (excluding orbits of baUistic reso- 
nances in the stochastic layer). All the trajectories in- 
side the CTB are ballistic, some of them are regular and 
the other ones are chaotic. The amplitude-frequency di- 
agram in Fig. [2] demonstrates clearly destruction of the 
CTB at some values of the perturbation parameters e 
and V and onset of cross-jet transport. 

Our Hamiltonian fiow with the streamfunction ^ is 
degenerate, i. e. it violates the non-degeneracy condition, 
df /dl ^ 0, for some values of the natural frequency of 



passive particles / and their actions / in the unperturbed 
system. Physically it means that the zonal velocity pro- 
file has a maximum. In theory of nontwist maps 
the curve, for which the twist condition (analogue of the 
non-degeneracy condition) is violated, is called a non- 
monotonic curve fH]. In our model fiow ^ it is some 
value of unperturbed streamfunction along which the fre- 
quency / is maximal (see Fig.[T](b)). 

Instead of integrating the advection equations |(T|) , we 
integrate the corresponding Poincare map, an orbit of 
which is defined as a set of points {{xi,yi)}l^_^ on the 
phase plane such that GT{xi,yi) = {xi+i,yi+i), where 
Gt is an evolution operator on a time interval t and T = 
2-1: 1 V is the period of perturbation. Operator Gt can be 
factorized as a product of two involutions Gt — hlo, 
where /i = Gt-^o is also a time reversal simmetry |37l]. 

A periodic orbit of period nT (n = 1, 2, . . . ) is an or- 
bit such that (xi+n, j/i+„) = (x^ + 27rm, y^), Vi, where 
TO is an integer. An invariant curve is a curve invari- 
ant under the map. The nonmonotonic curve is not an 
invariant curve under a perturbation. The winding (or 
rotation) number w of an orbit is defined as the limit 
w — lim [{xi — xo)/(27ri)], when it exists. The winding 

i — *C30 

number is a ratio between the frequency of perturbation 
V and the natural frequency /. Periodic orbits have ratio- 
nal winding numbers w = m/n. It simply means that a 
baUistic passive particle in the fiow files to frames before 
returning to its initial position xq (modulo 27r) after n pe- 
riods of perturbation. Winding numbers of quasiperiodic 
orbits are irrational. 

Now we are ready to introduce the important notion 
of a central invariant curve (CIC). We define CIC as an 
curve which is invariant under the operators S and Gt- 
It can be shown, that two curves invariant under S have 
at least two common points. The curves, which are in- 
variant under Gt, cannot intersect each other. So, the 
CIC is a unique curve. Following to Ref. [sHi, one can 
show that the CIC corresponds to a local extremum on 
the winding number profile with an irrational value of 
w. Such curves are called «shearless curves* in theory of 
nontwist maps [s^l • The significance of a shearless curve 
is that it acts as a barrier to global transport in the phase 
space of a nontwist map. The violation of the twist con- 
dition leads to existence of more than one orbit with the 
same winding number arising in pairs on both sides of 
the shearless curve. Those pairs of orbits can coUide and 
annihilate at certain parameter values. The collision of 
the orbits involves in phenomenon, which was called as 
reconnection of invariant manifolds of the corresponding 
hyperbolic orbits [3l|. 

The CIC should not be thought as the last cross-jet 
barrier curve in the CTB in the sense that it breaks 
down under increasing the perturbation amplitude in the 
last turn. Sometimes it is the case, but sometimes it is 
not. Nevertheless, the CIC serves a good indicator of the 
strength of the CTB and its topology. 

The CIC can be constructed by successive iterations 
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of so-called indicator points [s^. In our model flow ([2]) 
with the symmetries (0]) and JH) , indicator points are the 
points {x'p , Vj'^^)' ^ = li2, which are solutions of the 
equations 

or 

/,(xf,,f) = 5(.f,yf), (7) 

where index j numerates the points. The equation ^ 
gives a pair of indicator points: {x''^^ = tt/2, y[^'' — 0) 

and {x'i^ = 37r/2, y!^^ = 0). Instead of solving Eq. 
we solve the equivalent equation 

GT{x,y) ^ iQS{x,y) = {tt ~ x,~y). (8) 

If some {x, y) is a solution of ([8]), then Iq{x^ y) is a 
solution of The equation ^ cannot be solved ana- 
litically, so we apply the numerical method based on com- 
puting a minimum of the function r{x,y) = \ \GT{x,y) — 
{tt — x,—y)\\, where 1 1 • 1 1 is a norm on the cyHnder. Since 
y) > for any (a:, y), the points with r{x, y) — are 
minima of the function r{x, y). Thus, solution of Eq. ([8]) 
reduces to searching for a local minimum of the function 
r{x,y) with the additional condition r{x,y) — at the 
point of the minimum. There are a number of numer- 
ical methods for doing that job. We prefer to use the 
downhill simplex method. In our problem the function 
r{x, y) always has two minima, i. e. two indicator points 
transforming to each other under action of the operator 
S. 

Next, we study iterations, i. e. Poincare mapping 
of one of the indicator points (a;o,yo)- If the itera- 
tions {xi,yi) = &rp{xo,yo) are conflned between invari- 
ant curves in a bounded region, the following three cases 
are possible in dependence on the dimension d of the set 
{xi,yi): 

1) The iterations lie on a curve on the phase plane with 
d—1 which is a CIC. 

2) The iterations is an organized set of points with 
d = 0. It means that they constitute either a central 
periodic orbit or a central almost periodic orbit, an orbit 
that could not form a smooth curve on the phase plane 
for a limited integration time. 

3) The iterations form a central stochastic layer with 
rf = 2. 

If the iterations are not confined by any invariant 
curves in a bounded region, i. e. they occupy all the ac- 
cessible phase plane to the south and north from the cen- 
tral jet, then there exists global chaotic transport. Thus, 
the type of motion of indicator points provides an indica- 
tor of global chaos and the absence of barriers to cross-jet 
transport. 

Possible topologies of the corresponding CTB at the 
fixed frequency i/ = 1.2 and with increasing values of the 
perturbation amplitude e are illustrated in Fig. [3] plot- 
ting iterations of the indicator points computed by the 



above-mentioned method. The panel (a) illustrates the 
case when those iterations form a CIC. Another typical 
situation is shown in Fig. [3] (b) where the iterations fall 
in small segments filling at t ^ oo a continuous curve 
which is a central almost periodic orbit. If the iterations 
fill up not a curve but a bounded region between invari- 
ant curves, then there appears a central stochastic layer 
preventing cross-jet transport (Fig. [3] (c)). When itera- 
tions of the indicator points occupy a region that is not 
confined by any invariant curves, it means destruction 
of the CTB and onset of global chaos, i. e., chaos in a 
large region of the phase space accompanied by cross-jet 
transport. 

The indicator points have been found with the help 
of the above-mentioned numerical procedure, and their 
iterations have been computed in the following range of 
the control parameters: v € [0.95 : 1.5] and e € [0.01 : 1]. 
We assume that the iterations are bounded, if their co- 
ordinates do not cross the unperturbed separatrices after 
5 X 10"* iterations. The dimension d of the set of those it- 
erations is computed by the box-counting method, where 
the value of d for the box size = (1/2)'= is defined as 

dk^log,^, (9) 

where Nk is a number of boxes of the size ek containing 
set points. The dimension dk goes to zero with decreasing 
Cfc, and one cannot distinguish in this Hmit between the 
central almost periodic orbit and the central stochastic 
layer at large k. Comparing the values of dk at different 
values of k, we were able to find the empirical value k = A 
which is enough to make the difference. 

The results of computation of the dimension di{e,i') 
for a set of iterations of the indicator points are shown 
in the bird-wing diagram in Fig. [H That one and the 
other bird-wing diagrams in the parameter space show 
the properties of CTB and CIC in the range of compara- 
tively small values of the perturbation amplitude (0.01 < 
£ < 0.1) and the frequency (1.15 < v < 1.5) correspond- 
ing to particles moving in the central jet (Fig. [1] (b)). 
White color corresponds to the regime of global chaotic 
transport with unbounded motion of iterations of the in- 
dicator points. Otherwise, the CTB exists but its topol- 
ogy is different. Grey color means that there exists a 
CIC with 0.95 < ^4 < 1.05 in the corresponding range 
of the parameters. White rectangles, which are hardly 
visible in the main panel (see their magnification on the 
inset of the figure) , means existence of a central almost 
periodic orbit with ci4 < 0.95 and black strips — a central 
stochastic layer with c?4 > 1.05. 



B. Geometry of the central invariant curve and its 
bifurcations 

To quantify complexity of the CIC we define its length 
L as a sum of the distances between the iterations of the 
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indicator points (x^, y^) ordered on the phase plane in the 
following way: 

1. The first step. A point Bq, belonging to a set of 
iterations of the Poincare map {xi,yi), is marked. 

2. The (j+l)-th step. We find and mark among all the 
unmarked points that one, -B^+i, which minimizes 
the EucHdean distance Dj = D{Bj, Bj+i) between 
Bj and -B^+i. 

3. The procedure is repeated unless all the points will 
be marked. 





2 * ^ ' * 



FIG. 3: Poincare mapping of indicator points. In the 
first three panels the orbit of these points is bounded and 
there exists a central transport barrier (CTB) with (a) GIG 
(e = 0.01, — 1.2), (b) central almost periodic orbit (e = 
0.011997277, u = 1-2), and (c) central stochastic layer at 
(e = 0.01177, u — 1.2) with the inset demonstrating a magni- 
fication of a small region, (d) Destruction of CTB and onset 
of global chaotic transport as a result of unbounded iterations 
of indicator points (e = 0.041, f — 1.2). 

As an output we have an ordered set of points Bj 
constituting a CIC. The accuracy is controlled by the 
quantity maxDj. Large values of this quantity mean 
that the points are ordered in a wrong way or a set of 
points is chaotic. To increase the number of points we 
use in addition to the original points {xi,yi) their "im- 
ages" {xi +Tr, —yi) as well. To minimize the computation 
time the points are sorted in accordance with their x co- 
ordinates. 

Figure [5] illustrates metamorphosis of the CIC as the 
perturbation amplitude increases. We start with the 
CIC, shown in Fig. [5] (a), which we call a nonmeandering 
CIC. At the critical value e ~ 0.011758, invariant man- 
ifolds of hyperbolic orbits of two chains of the 1:1 reso- 
nance islands on both sides of the CIC connect, and after 
that the CIC becomes a meandering curve of the first or- 
der (Fig.[5](b)) and period T. The period of CIC's mean- 
dering is simply a period of nearby main islands [35l.[4l]|. 
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FIG. 4: Bird-wing diagram of the box-counting dimension 
d4,{e,v). White color: regime with global chaotic transport 
with unbounded motion of indicator points (see Fig. [3] (d)). 
If the motion of indicator points is bounded, then there exists 
a GTB but its topology may differ. Grey color (0.95 < < 
1.05): regime with a CIC (see Fig.[3](a)). Small white regions 
which are hardly visible inside the grey "wing" {d^ < 0.95): 
regime with a central almost periodic orbit (see Fig. [3](b)). 
Black color (d4 > 1.05): regime with a central stochastic layer 
(see Fig. [3] (c)). Inset shows magnification of a small region 
in the parameter space with visible white and black regions. 



At the next critical value e — 0.01178721, reconnection of 
invariant manifolds of secondary resonance islands takes 
place. The corresponding second-order meandering CIC 
with period 79T is shown in Fig. [5](c). Highly meander- 
ing CICs of higher orders appear with further increasing 
the perturbation amplitude (Fig. [5] (d)). 

Some smooth invariant curves inside the CTB break 
down under the perturbation ^ , and chains of baUistic 
resonance islands appear at their place. Those islands 
appear in pairs to the north and south from a CIC due to 
the fiow symmetries ^ and ^ (see Fig. [5] (a) and (b)). 
Geometry of the CIC, size and number of the islands, 
and topology of their invariant manifolds change with 
variation of the perturbation amplitude e and frequency 
J' in a very complicated way. 

In Fig. [6] we plot in the parameter space the values of 
the CIC length L coding it by nuances of the grey color. 
White color corresponds to those values of the parame- 
ters £ and ly for which cross-jet transport exists due to de- 
struction of the CTB. Black color codes the regime with a 
broken CIC but a remaining CTB that prevents cross-jet 
transport (^4 > 1.05). Dotted and dashed lines on the 
plot are the resonant bifurcation curves along which the 
CIC winding number w is rational. The m/n resonant 
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FIG. 5: Metamorphosis of the central invariant curve (CIC). 
(a) Nonmeandering CIC {u = 1.2, e = 0.01174929). (b) 
Meandering CIC of the first order and period T (v = 1.2, 
e = 0.01178721). (c) Meandering CIC of the second order 
and period 79T {u = 1.2, e = 0.01179027). (d) Meandering 
CIC of a higher order {u = 1.2, e = 0.01179339). 



bifurcation curve is the set of values of the control pa- 
rameters for which a reconnection of invariant manifolds 
of the n : m resonances takes place. The dotted lines cor- 
respond to even resonances with w = {2k — 1) / 2fc and 
the dashed lines are odd resonances with w = 2kl (2fc-|-l), 
k = 1,2,.... All those curves end up in the dips of the 
bird- wing diagram. 

In order to analyze a fractal-like boundary of the bird- 
wing diagram in Fig. [6l we cross it horizontally at the 
frequency v = \.2 and consider the plot in the range 
of interest of e (Fig. [7] (a)). The perturbation frequency 
1/ = 1.2 is close to the maximal frequency /max of parti- 
cles in the middle of the jet in the unperturbed flow (see 
Fig. [1] (b)). The plot L(e) consists of a number of spikes 
with different height and width. In the range of small 
values of the perturbation amplitude (e < 0.011756), the 
length of the CIC is approximately the same L k, 7.35 (a 
small fragment of the function L{e) is shown in Fig. [7] (a) 
just to the left from the vertical line 1/1). In that range, 
the CIC is a nonmeandering curve (see Fig. [5] (a)) sur- 
rounded by smooth invariant curves and 1 : 1 resonance 
islands with a heteroclinic topology. The size of those 
islands is comparable with the frame size. The width 
of the CTB, filled by invariant curves around the CIC, 
decreases with increasing the perturbation amplitude e. 

The CIC winding number w changes under a varia- 
tion of the perturbation amplitude. At e ~ 0.011756, 
invariant manifolds of the 1:1 resonance connect, and a 
central stochastic layer appears at the place of the CIC. 
This layer exists up to e ~ 0.011785 (see a random set of 
points in Fig. [7| (a) in that range of e). At e > 0.011785, 
the CIC appears again. Now it is a meandering curve 
of the first order (see Fig. [5] (b)) whose length is larger 
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FIG. 6: Bird-wing diagram showing the length L of the CIC 
by nuances of the grey color in the parameter space {e,v). 
White zone: regime with a broken CTB and cross-jet trans- 
port. Black color: regime with a broken CIC but a remain- 
ing CTB preventing cross-jet transport. Resonant bifurcation 
curves, along which the CIC winding numbers w are rational, 
end up in the dips of "the wing". Dotted and dashed lines 
correspond to even and odd resonances, respectively. 

due to reconnection of the 1:1 resonance islands. As e 
increases further, the CIC length L changes in a wide 
range. Smooth fragments with approximately the same 
value of L « 20 alternate with spikes of different height 
and width. The spikes are condensed, when approach- 
ing to the value w — 1/1, and overlap in the range 
e ~ [0.011756 : 0.011785]. 

The arrangement of the spikes in Fig. [7] (a) can be 
explained using a representation of rational numbers by 
continued fractions. A continued fraction is the expres- 
sion 

c = [ao;ai,a2,a3, . . .] = ao H , (10) 

ai H J 

a2 H 

as + • • • 

where oq is an integer number and the other a„ are nat- 
ural numbers. Any rational (irrational) number can be 
represented by a continued fraction with a finite (infi- 
nite) number of elements. The spikes in Fig. [7] (a) are 
arranged in convergent series in such a way that each 
spike in a series generates a series of spikes of the next 
order. For example, the series of the integer n : 1 reso- 
nance has the winding numbers equal to 1/n or [0;ri] in 



0.01193 ^sc 0.011935 g 

FIG. 7: Dependence of the length of GIG L on the pertur- 
bation amplitude at the fixed frequency v — 1.2. (a) Gen- 
eral view of the dependence L{e) in the range of interest 
£ — [0.01175 : 0.012]. Vertical dotted lines correspond to 
rational values of the GIG winding number w. The resonance 
1/1 appears at £ = 0.11756. The arrangement of the spikes 
is explained in the text, (b) Magnification of one of the wide 
spikes in panel (a). Solid line is a winding number profile w{e) 
with the value w = 62/63 shown by the dashed line. 




FIG. 8: Poincare sections at f = 1.2 and increasing values of £ 
in the range corresponding to the wide spike with w = 62/63 
in Fig. [7](b). (a) Meandering GIG surrounded by meander- 
ing invariant curves {e = 0.011931). (b) GIG meandering 
between the odd 63 : 62 islands born as a result of a saddle- 
center bifurcation (£sc = 0.011934). (c) GIG destruction due 
to connection of invariant manifolds of the 63 : 62 islands 
(£ = 0.01193511). A narrow stochastic layer appears at the 
place of the GIG. (d) GIG appears again {e = 0.0119352). 



FIG. 9: Poincare sections at £ = 0.015 and increasing val- 
ues of u. (a) GIG between islands of the even 2:1 resonance 
{u = 2.51 < 2/inax = 2.556). (b) Reconnection of invari- 
ant manifolds of that resonance and a formation of a vor- 
tex pair with a narrow stochastic layer shown by bold curves 
{u = 2.55). (c) The vortex size decreases with increasing v 
{u = 2.555). (d) Past some critical value of the vortex pair 
disappears and GIG appears again {v — 2.56). 



the continued-fraction representation. Each spike in that 
series generates a series of resonance spikes of the next 
order converging to the parent spike. Winding numbers 
of those resonances are [0; n, i], i = 2, 3, 4 . . . (at i = 1, 
one gets a spike in the main series because of the iden- 
tity [ao;ai, . . . , a„, 1] = [ao; ai, . . . ,a„ -I- 1]). The spikes 
with [0; 1, i] = i/(i -I- 1) converge to the spike of the 1 : 1 
resonance. That is clearly seen in Fig. [7| (a). The di- 
rection of convergence of the spikes in a series alternate 
with the series order: the winding number increases in 
the series of the first order, decreases in the series of the 
second order, and increases again in the series of the third 
order. That is why a chaotic region in Fig. [7] (a) is sit- 
uated to the right from the 1 : 1 resonance, i. e. in the 
range of smaller values of w. Whereas, it is to the left 
for the series of the second order, i. e. in the range of 
larger values of w. There also exists an additional hier- 
archical structure with fractional 1 : n resonances, whose 
frequencies are below the 1 : 1 resonance frequency, and a 
series with resonances corresponding to the spikes below 
[0; 1, i, (1)], for example, a clearly visible series of spikes 
below [0; 1,4,(1)] converging to the spike 5/6 in the L 
diagram (see Fig. [6]). 

Unfortunately, we could not identify series of the third 
and a higher order because numerical errors in identify- 
ing the winding numbers are greater than the distance 
between the spikes of higher-order series. Unresolved re- 
gions on the plot L(e) in Fig. [7] (a) appear because of a 
decrease of the distance between the spikes in the same 
series with increasing series number, a process resembhng 
Chirikov's overlapping of resonances. 
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A magnification of one of the wide spikes is shown in 
Fig. El (b). We plot the winding number profile w{e) to- 
gether with the function L(e) for the spike. To illustrate 
what happens with the CIC and its surrounding with in- 
creasing £ we plot the corresponding Poincare section in 
Fig. [a In the range e ~ [0.0119 : 0.01192] the lengths of 
the CIC and surrounding invariant curves increase slowly 
due to small changes in their geometry (Fig. [8] (a)). Af- 
ter a saddle-center bifurcation at Egc — 0.011934, there 
appear two chains of homoclinically connected 63 : 62 
islands separated by a meandering CIC (Fig.[8](b)). The 
amplitude of the CIC meanders increases with further 
increasing e in the range e ~ [0.011934 : 0.011935]. 
In that range the CIC disappears and appears again in 
a random-like manner (see the corresponding fragment 
on the plot L(e) in Fig. [7] (b)) due to overlapping of 
higher-order resonances and reconnection of their invari- 
ant manifolds. The example of such a reconnection for 
the 63 : 62 resonance at e = 0.01193511 is shown in 
Fig. [8] (c) where a stochastic layer appears at the place of 
the CIC. As £ increases further, the CIC appears again 
but with a smaller number of meanders (Fig. [8] (d)). 
Animation of the corresponding patterns is available at 
[http: / / dynalab.poi.dvo.ru / papers /c ic.avij 

The other wide spikes in the plot L{e) with a similar 
structure are caused by another odd resonances between 
the external perturbation and particle's motion along the 
CIC. Under a CIC resonance with the winding number 
w = m/n,we mean reconnection of invariant manifolds of 
the resonance n : m and onset of a local stochastic layer. 
The narrow spikes, situated between the wide ones in 
Fig. El (a) , correspond to reconnection of even resonances. 
They are hardly resolved on the plot. Even resonances 
of higher orders have a smaller efifect on CIC geometry 
then odd resonances. As an example, we illustrate in 
Fig. [9| metamorphosis of the CIC with the winding num- 
ber w = 1/2. The perturbation amplitude is fixed at a 
rather small value £ — 0.015 and the frequency increases 
in the range 2.51 < < 2.556 = 2/max. At ^ = 2.51 
there are islands of the even 2 : 1 resonance separated by 
a CIC (Fig. m (a)). At some critical value of v invariant 
manifolds of the 2 : 1 resonance connect and the islands 
form a tight vortex-pair structure surrounded by a nar- 
row stochastic layer (see Fig.[9|(b) at v = 2.55). The size 
of the pair decreases gradually with further increasing v, 
and the corresponding hyperboHc orbits approach each 
other (see Fig.[9|(c) at i' = 2.555). At some critical value 
of hyperbolic and elliptic orbits of the resonance col- 
lide and annihilate, and CIC appears again (see Fig.[9|(d) 
at u = 2.56). Vortex pairs of the other even resonances 
are formed in a similar way. The higher is the order of 
the resonance, the smaller is the vortex size. 

We conclude this section by computing winding num- 
bers w of the CIC in the parameter space. The result is 
shown in the bird- wing diagram in Fig.[TOl The CIC does 
not exist in the white region where the CTB is broken 
and cross-jet transport takes place. The curves, which 
end up on the tips of the "feathers of the wing", have 



winding numbers w with the following continued-fraction 
representation: [ao; ai, . . . , o„, (1)]. These are so-called 
noble numbers which are known to be the numbers that 
cannot be approximated by continued-fraction sequences 
to better accuracy than the so-called Diophantine condi- 
tion (see, for example, [i2|). The CICs with noble wind- 
ing numbers are in a sense the most structurally robust 
invariant curves, i. e. they may survive under a compar- 
atively large perturbation preventing cross-jet transport. 
The noble curves are arranged in series like the resonant 
bifurcation curves with rational winding numbers which 
end up in the dips of "the wing" in the bird-wing dia- 
gram in Fig. m For example, the noble series [0; 1, z, (1)] 
in Fig. [ini corresponds to the resonance series [0;l,i] in 
Fig. m In the w diagram we show a few representatives 
of the noble series [0; 1, i, (1)] and series of the next order 
[0; 1, z, j, (1)] (see Fig. M with j = 2, 3). 
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FIG. 10: Bird-wing diagram w{e,u) in the parameter space 
showing values of the winding number w of the CIC by nu- 
ances of the grey color. White zone: regime with broken CTB 
and cross-jet transport. The curves with irrational winding 
numbers end up on the tips of the "feathers of the wing" (some 
of them are marked by the corresponding noble numbers), 
whereas the curves with rational winding numbers (shown in 
Fig. [6]) end up in the dips of "the wing". 



IV. BREAKDOWN OF CENTRAL TRANSPORT 
BARRIER 

We have studied in the preceding section properties of 
the CIC which has been shown to be a diagnostic means 
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FIG. 11: Destruction of central transport barrier upon moving 
in the parameter space along a resonant bifurcation curve with 
the rational winding number u; = 8/9. (a) Narrow stochastic 
layer on the Poincare section is confined between invariant 
curves which provide a transport barrier. The perturbation 
parameters (e = 0.04889, ly' = 1.31625) are chosen on the 
curve with w = 8/9 nearby its right edge (see Fig. [GJ. (b) 
Onset of cross-jet transport at the values of parameters (e — 
0.054, u = 1.285) chosen in the white zone of that dip. 
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FIG. 12: Destruction of central transport barrier upon mov- 
ing in the parameter space along a curve with the noble value 
of the CIC's winding number w — [0; 1,5, 1, (1)]. When ap- 
proaching a tip of the corresponding "feather of the wing" 
in Fig. IIOI one observes on the Poincare section a decrease 
in the width of transport barrier with a CIC (bold curves) 
inside, (a) e = 0.07017, v = 1.367875. (b) e = 0.07416, 
V = 1.350375. (c) e = 0.0796067, ly = 1.325875. (d) Onset 
of cross-jet transport at the values of parameters (e = 0.08, 
v = 1.3142) chosen beyond a tip of that "feather" in the white 
zone in Fig. 1101 



to characterize CTB and its destruction. CTB separates 
water masses to the south and the north from the cen- 
tral jet and prevents their mixing. It is not a homoge- 
neous jet-Uke layer but consists of chains of ballistic is- 
lands, narrow stochastic layers, and meandering invariant 
curves of different orders and periods (including a CIC) 
to be confined by invariant curves from the south and 
the north. Those curves break down one after another 
when increasing the perturbation amplitude e, produc- 
ing stochastic layers at their place on both sides of the 
central jet, until the stochastic layers merge with one an- 
other and with stochastic layers around the southern and 
northern circulation cells producing a global stochastic 
layer and onset of cross-jet transport. 

Upon moving along any resonant bifurcation curve 
with a rational value of the winding number w in the 
bird-wing diagram in Fig. [6l we have those values of the 
perturbation amplitude e and frequency ly at which the 
corresponding CIC is broken due to reconnection of in- 
variant manifolds. It does not mean that CTB is bro- 
ken as well. That is the case only if we are at the dips 
of the "wing". The process of CTB destruction for this 
type of movement in the parameter space is illustrated 
in Fig. El We fix a point (e = 0.04889, ly = 1.31625) on 
the resonant bifurcation curve with w = 8/9 nearby its 
right edge in Fig. [6] and plot the corresponding Poincare 
section. A narrow stochastic layer, confined between in- 
variant curves providing a transport barrier, appears on 
the Poincare section in panel (a) at the place of a broken 
CIC. The barrier will be broken if one would choose the 
values of parameters in the white zone in Fig.[6l Merging 
of southern and northern stochastic layers and onset of 
cross-jet transport are shown in panel (b) at e = 0.054, 
ly = 1.285. 

Upon moving along any curve with a noble value of the 
winding number w in the bird-wing diagram in Fig. [TOl 
we have those values of the perturbation amplitude e and 
frequency ly at which a CIC with the corresponding noble 
number exists. The process of CTB destruction for the 
motion in the parameter space along the noble curve with 
w = [0; 1, 5, 1, (1)] is illustrated in Fig.[l2l When moving 
to the tip of the corresponding "feather of the wing" in 
Fig. [ini one observes progressive destruction of invariant 
curves and decrease of the width of the transport barrier 
(panels (a) and (b)) unless a single CIC remains as the 
last barrier to cross-jet transport (panel (c)). Onset of 
cross-jet transport (panel (d)) happens at the values of 
parameters chosen beyond a tip of that "feather" in the 
white zone in FigfTOl 

Upon moving along any resonant bifurcation curve to 
the corresponding dip of the bird- wing diagram in Fig.[6l 
we find cross-jet transport at smaller values of the pertur- 
bation amplitude as compared to the case with irrational 
winding numbers because in order to provide cross-jet 
transport in the first case it is enough to destruct all the 
KAM curves. Whereas, CICs with irrational and espe- 
cially noble values of the winding number may deform in 
a complicated way but still survive under increasing e up 
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to comparatively large values. 

V. CONCLUSION 

Being motivated by the problem of cross-jet transport 
in geophysical flows in the ocean and atmosphere, we 
have studied in detail topology of a central transport 
barrier (CTB) and its destruction in a simple kinematic 
model of a meandering current with chaotic advection of 
passive particles (Fig. [TJ that belong to the class of non- 
degeneracy Hamiltonian systems. Direct computation of 
the amplitude-frequency diagram (Fig. [2|) demonstrated 
onset of cross-jet transport at surprisingly small values of 
the perturbation amplitude e provided that the perturba- 
tion frequency v was sufficiently large. As an indicator 
of the strength of the CTB and its topology, we used 
a central invariant curve (CIC) which was constructed 
by iterating indicator points by a numerical procedure 
borrowed from theory of nontwist maps. The CTB has 
been shown to exist provided a set of the iterations was 
bounded. Otherwise, cross-jet transport has been ob- 
served (Fig. [31). The results were presented as a diagram 
of the box-counting dimension of those sets of iterations 
in Fig. H 

Geometry of the CIC has been shown to be highly sen- 
sitive to small variations in the parameters near a fractal- 
like boundary of the diagram (Fig. [5]). Quantifying com- 
plexity of the CIC's form by its length L, we computed 
the corresponding L diagram looking Hke a bird wing 
with a fractal-like boundary (Fig. [6]). Resonant bifurca- 
tion curves with rational winding numbers mjn end up in 
the dips of the boundary. Along those curves in the pa- 
rameter space, invariant manifolds of the corresponding 
n : m resonances connect providing a destruction of the 
CIC. Scenarios of the reconnection are different for odd 
(Fig. [HI) and even (Fig. [9]) resonances. Animation of the 
process at http:/ /dynalab. poi.dvo.ru/papers/cic.avi pro- 
vides a visual demonstration of complexity of topology 
of the CTB and its destruction. Computing the winding 
number w of the CIC, we have got an information about 
those values of the perturbation parameters at which the 
CTB is strong or weak. Using representation of the val- 
ues of winding numbers by continued fractions, we were 
able to order spikes with rational values of w in Fig.[7|into 
hierarchical series of the corresponding CIC resonances. 
The curve, which end up on the tips of "feathers of the 
wing" in the winding-number diagram in Fig. [TOl have 
noble winding numbers which are so irrational that the 
corresponding CICs break down in the last turn when 
varying the perturbation parameters. The noble curves 
have been found to be arranged in series like the resonant 
bifurcation curves with rational values of w. Destruction 
of CTB is illustrated for two ways in the parameter space: 
upon moving along resonant bifurcation curves with ra- 



tional values of w (Fig. \TT\\ and along curves with noble 
values of w (Fig. [T2|) . 

In conclusion we address two points that may be im- 
portant in possible aplications of the results abtained. 
Molecular diffusion in laboratory experiment and turbu- 
lent diffusion in geophysical flows are expected to wash 
out ideal fractal-Hke structures caused by chaotic advec- 
tion after a characteristic time scale. The question is 
what is this scale. As to molecular diffusion in the ocean, 
the diffusion time-scale Li^ /D is very large since the dif- 
fusion coefficient is of the order of D ~ 10~^ cm^/sec and 
L is of a kilometer scale. The scale of molecular diffusion 
in laboratory tanks is, of course, much smaller. How- 
ever, some fractal-Hke structures have been observed in 
real laboratory experiments (see, for example, Ref. ji^l 
and the book |[1^] for a recent review of experiments). 
Modelling of a combined effect of chaotic advection and 
turbulent diffusion in the ocean is a hard problem de- 
serving a special consideration. Any kind of diffusion is 
expected to intensify cross-jet transport. 

The advantage of the kinematic approach is its ability 
to identify different factors that may enhance or sup- 
press cross-jet transport. However, the results obtained 
with our simplified kinematic model should be taken with 
caution to describe cross-jet transport in real geophysical 
flows. In any kinematic model the velocity fleld is postu- 
lated based on known features of the current while in dy- 
namic models it should obey dynamical equations follow- 
ing from the conservation of potential vorticity (il. [45|. 
It is very difficult to formulate an analytic and dynami- 
cally consistent model with chaotic advectio n ( for a dis- 
cussion and examples of such models see [Tsl. [27l. [29l. li^l) . 
One approach is to seek solutions of the fiuid dynam- 
ics equations that are self-consistent to linear order [47| . 
Some aspects of cross-jet transport in a linearized model 
with a zonal Bickley jet current and two Rossby waves 
have been studied in Refs. [13, HI, US] • Potential vorticity 
is not exactly conserved within the linear approximation 
but models that are self-consistent to linear order provide 
a compromise between the self-consistency demands and 
the fruitfulness of Hamiltonian models. 

The question, how predictions of kinematic models in 
destruction of barriers to cross-jet transport carry over 
to more realistic dynamical models, remains open. We 
plan in the future to apply the methods developed in 
the present paper to a dynamically consistent model of a 
meandering current with Rossby waves. 
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